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tJ> Abstract 

This paper is devoted to study the effects of electromagnetic field 
on the energy density inhomogeneity in the relativistic self-gravitating 
fluids for spherically symmetric spacetime. Two important equations 
of the Weyl tensor are formulated which help to analyze the energy 
density inhomogeneity in this scenario. We investigate two types of 
fluids, i.e., non-dissipative and dissipative. The non-dissipative fluid 
further includes dust, locally isotropic, and locally anisotropic charged 
fluids. We explore the effects of different factors on energy density 
inhomogeneity in all these cases, in particular, the effect of charge. 
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1 Introduction 

Gravitational collapse emerges from a highly inhomogeneous initial state 
which can be explained in terms of inhomogeneous energy density distri- 
bution. Energy density inhomogeneity plays a vital role in the collapse of 
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self-gravitating fluid. Penrose [I] discussed the role of the Weyl tensor in 
the evolution of self-gravitating system. He provided a simple relation be- 
tween the Weyl tensor and energy density to investigate the gravitational 
arrow of time. The Weyl tensor may be represented exclusively in terms 
of energy density and local anisotropy of pressure, which affects the fate of 
gravitational collapse. 

There is a large body of literature available [2]-|E] which indicate the im- 
portance of energy density inhomogeneity in self-gravitating fluid. Joshi and 
Dwivedi [2] analyzed the Tolman-Bondi model to examine the nature and 
occurence of naked singularities for inhomogeneous gravitational collapse. 
Triginer and Pavon [3] discussed heat transport equation in an inhomoge- 
neous spherically symmetric universe. Mena and Tavokol [I] investigated 
how it is important in the dust collapse. Herrera et al. j5] derived a relation 
for the active gravitational mass of collapsing fluid distribution to explain the 
effects of density inhomogeneity and local anisotropy on spherical collapse. 
The same authors [6J also studied the behavior of locally anisotropic, self- 
gravitating spherical symmetric dissipative fluid with the Weyl tensor and 
density inhomogeneity. 

The study of self-gravitating spherically symmetric charged fluid distribu- 
tion has been the subject of interest for many people. Rosseland [7] and Ed- 
dington [8] started the study of self-gravitating spherically symmetric charged 
fluid distribution. Since then, a lot of work has been done to investigate the 
effects of electric charge on the structure and evolution of self-gravitating 
system [9]-[T6]. Di Prisco et al. [17] investigated the non-adiabatic charged 
spherically symmetric gravitational collapse as well as the energy density in- 
homogeneity. Sharif and his collaborators [18] derived dynamical as well as 
transport equations of matter dissipating in the form of shear viscosity to see 
the effect of charge on gravitational collapse. 

The gravitational collapse is highly dissipative process [HI [20], hence its 
effects are much important in the study of collapse. Misner [21] explained 
the dissipative relativistic gravitational collapse by using the streaming out 
approximation. Two possible approximations (diffusion and streaming out) 
are usually considered in dissipative process. Diffusion approximation is the 
approximation for which energy flux of radiation (like thermal conduction) 
is proportional to the gradient of temperature. Israel and Stewart [22J for- 
mulated transport equation required for diffusion approximation . Lattimer 
[23] explained that during emission process, the role of radiation transport is 
closer to diffusion approximation than the streaming out approximation. Re- 
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cently, Herrera [24J explored energy density inhomogeneity in self-gravitating 
fluid as well as its stability. 

In this paper, we extend this work to see the effects of charge on en- 
ergy density inhomogeneity due to fluid distribution. We shall discuss the 
effects of electric charge on the relationship between the Weyl tensor and 
energy density inhomogeneity. The outline of the paper is as follows. In 
the next section, we review the relevant kinematics and Einstein-Maxwell 
field equations for a spherically symmetric distribution of collapsing charged 
fluid. Section 3 is devoted to dynamical as well as transport equations and 
two equations of the Weyl tensor. In section 4, different aspects of fluid 
are considered to explore energy density inhomogeneity. The last section 
contains the conclusion of the results. 



2 Spacetime and Matter Distribution 

We take spherically symmetric distribution of dissipative collapsing charged 
fluid bounded by a spherical surface S in the comoving coordinates as 

ds 2 = -A 2 dt 2 + B 2 dr 2 + C 2 (d6 2 + sin 2 (1) 

where A, B and C are functions of t and r. Matter under consideration is 
anisotropic fluid suffering dissipation in the form of heat, i.e., 

= (A* + P±)V a V p + P±g aP + (P r - P±)x«Xp + q a Vp + V a q p + djp, (2) 

where /i, P r , P±, V a , x a , q a , l a and e represent the energy density, the radial 
pressure, the tangential pressure, the four velocity of the fluid, a unit four- 
vector along the radial direction, the heat flux, a radial null four-vector and 
the energy density of null fluid describing dissipation in the free streaming 
approximation respectively. These quantities satisfy 

V a Vp = -l, V a qp = 0, X°X/3 = 1, i a v a = -i, x a v a = o, l a l a = 0. 
Also, we define 

V a = A- 1 5%, q a = qB- 1 5^, r = A- 1 6% + B- 1 6?, x a = B~ x 5^. 
The shear tensor is given by 
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where h a p = g a p + V a Vp. The four acceleration a a and the expansion are 

a a = V a , p V^ Q = V a ;Q , 

which yield 

The non-zero components of the shear tensor will be 

2 D 2 °33 1 ^2 

Cll = °"> °"22 = . 2a = ^ 

o sin 6/ o 

with its scalar 



The electromagnetic energy-momentum tensor is 

T^ m) = ^(^7 - (5) 

where F a( g is the electromagnetic field tensor. In terms of four- vector formu- 
lation, the Maxwell field equations can be written in the form 

F a p = <f>p, a - = VoJ a , (6) 

where <f> a , J a and p represent the four potential, four current and magnetic 
permeability, respectively. Since charge is considered at rest in comoving 
coordinates, so J a and <p a become 

K = r = P v a . 

Here the charge density p and electric scalar potential 4> are both functions 
of t and r. The charge conservation equation gives 

rr 

s( r ) = 4tt / pBC 2 dr, (7) 
Jo 
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where s(r) is the electric charge interior to radius r. For metric (0Q), the 
Maxwell equations turn out to be 
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Equation (jSJ) yields 
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which satisfies Eq.(j9]). The field equations for charged dissipative fluid, G a p 
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3 Dynamical and Transport Equations 

Using the Misner-Sharp definition [25], the mass function in the presence of 
charge is given by 



CP 
2 



m = ^R 23 23 + — = ^ 



C 

2C ~ 2 



C 



h (15) 



The proper time and radial derivatives are defined by [T7] 



Id Id 



(16) 



where R stands for proper areal radius. Using this definition, the velocity of 
the collapsing fluid is U = D T C = S < 0. Thus we can write Eq. (fT5]) as 



;i7) 



Define jl — ji + e, P r = P r + e, q = q + e. The rate of change of mass is 



Dtttl = —4ir 



P r U + qE 



C 2 , D R m = 4tt 



fi + q 



E 



C 2 + -D R s. (18) 



This shows how mass is affected by different quantities. Integration of the 
second equation yields 
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The contracted Bianchi identities, ^( m ) Q,/3 4- J"( 



em) 



a/3 N 



0, yield 



_|_ j.( 
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(20) 



0, (21) 



where U = P r — P±. These are called dynamical equations. 
The Weyl tensor is defined by 
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The electric part of the Weyl tensor, E a/3 = C a ^ u V^V u , gives 



2 1 

E\\ = -B e, E 2 2 = — -C e, E : 
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while its magnetic part vanishes due to spherical symmetry. We may also 
write E a p as 

E a p = s(xaXf3 ~ -ha?)- (24) 
Solving Eqs.([IID-([l3D with Eqs. (|I5jl and ([23]), we get 



3m 2s 2 
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Using Eqs. flTS]) and (f2"3"|) . we can write 

C 

[e - 4tt(/} - n)] = 3- 47r(/i + P ± ) 
[e - 47r(/i - IT)]' = -3— ; Ivrll + 



c 



2s 2 
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(26) 
(27) 



These equations yield a relationship between the Weyl tensor, energy density 
and charge that help to discuss energy density inhomogeneities given in the 
next section. 

The corresponding transport equation for heat flux derived from Muller- 
Israel- Stewart theory [T8], [26] is given as 



1 frVP\ 

Th^\nq p . a + q a = -Kh a \T )fi + Tap) - -kT 2 ^— J q a , (28 



where k, T and r represent thermal conductivity, temperature and relax- 
ation time, respectively. Due to symmetry, the above equation has only one 
independent component given by 



rq 



T qT2 



l -rqQA-^{TA)'-qA. 



(29) 



For r = 0, Eckart-Landau [27] equation is recovered. In case of truncated 
version of the theory, the last term in Eq. fl28|) is removed and finally, we get 



rq 



K 

B 



(TA)> - qA, 



(30) 



which is the same as given in [21] . This shows that the electromagnetic field 
does not affect the heat transport equation. 



4 Energy Density Inhomogeneity 

In this section, we consider different aspects of fluid distribution that are 
responsible for energy density inhomogeneity . 
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4.1 Non-dissipative Charged Dust 

Firstly, the case of non-dissipative charged dust is taken, i.e., q = P r = P±_ = 
e = 0. Since the fluid is moving along geodesic (i.e., a a = = a), so Eq.flSJ) 
leads to A — 1. Consequently, Eqs. fl26]) and f l27|) reduce to 



(e - AM + 3^ (e - 4^ + ^\ = 0, (31) 
(e-An^ = -3^s-2^ + S ^. (32) 



For s = 0, this case reduces to the uncharged case [21]. Equation ( 1321) 
implies that for e = = s, we get // = 0. This shows that energy density 
inhomogeneity depends not only on the Weyl tensor but also on charge s. 
For e = = // (i.e., when spacetime is conformally flat and energy density 
vanishes), it follows from Eq.( l32l that C 2 (t,r) = s(r)0(t), where 0(i) is an 
arbitrary function. When // = 0, we obtain 

1 r,ss> 2s 2 C> 

-)dr, (33) 



C 3 7 v C C 2 

where integration function is chosen such that e(t, 0) = 0. This equation 
shows that homogeneity in energy density implies the existence of e in the 
presence of electromagnetic field. 

Next, we make use of Eqs.(jlj) and (120]) in (|3T|) so that 

Here integration function is chosen as e(0, r) = 0. With the help of Ray- 
chaudhuri equation and the field equations, we get an evolution equation for 
shear (see [?] for detail]) given by 

i+ °-l + 2 f = -e. (35) 

Equations ( 134"1) and ( 1351) show that conformal flatness and shearfree condi- 
tions do not imply each other. In this case, the shearfree implies vanishing 
of the Weyl tensor but converse is not true. 
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4.2 Locally Isotropic Non-dissipative Charged Fluid 

This case corresponds to a non-dissipative charged isotropic fluid, i.e., II 
q = € = 0, P r = P ± = P. Equations ([26} and (EZJ) reduce to 



(e - 4tt/*)- + 3^ - 4vr(/i + P) + = 0, (36) 

(s - 4^)' = -3^£ - 2^ + ^. (37) 

Equation (1371) is exactly the same as (1321) and hence the same behavior for 
e = q = 0, e = // = and // = 0. Using Eqs.gD and ((2DJ| in Q35), we have 

e = ~ y o [(^ + PMaC 3 ]^+^, (38) 

where e(0, r) = 0. For non-dissipative locally isotropic fluid (using Ray- 
chaudhuri equation as well field equations as in the non-dissipative charged 
case), the evolution equation of the shear takes the form 



a a 9 a 



2 



2„ C 



It is important to mention here that neither the vanishing of the Weyl tensor 
implies shearfree nor the shearfree condition corresponds to conformally flat- 
ness (i.e., the Weyl tensor disappears). If we assume the fluid to be shearfree, 
then Eq.(|38]) yields 

2s 2 

(40) 

This shows that under the effect of electromagnetic field for initial homoge- 
neous configuration, i.e, e(0,r) = 0, the Weyl tensor does not disappear at 
any time t. 

4.3 Locally Anisotropic Non-dissipative Charged Fluid 

In this case, we consider the role of pressure anisotropy with q = e = but 
II / 0. Under these conditions, Eqs.()26l) and ( 1271) become 

(e - 4tta* + 47rn)- + 3^ (e - 4vr(^ + P ± ) + ^\ = 0, (41) 
(e - 4^ + AnU)' + 3^ (e + ArrU + ^ - = 0. (42) 
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In the case of locally anisotropic fluid, the quantity e + Anil plays an im- 
portant role instead of the Weyl tensor. The existence of the energy density 
inhomogeneity depends on e + 47riT as well as charge unlike previous cases. 
From Eq. fj42|) . it follows that for e + AttU = = s, we get // = 0. When we 
take e + AirU = = //, then Eq.flU) reduces to C 2 (t, r) = s(r)0(t) while for 
/jf = 0, this yields 



e-MI + -j + 3 ^S— = - + —. (43) 

Using Eqs.g]) and (E01) in (HT)). it follows that 



-e - 47iTI + — + 3 



s 2 \ , o (-e-4vrn + ^)C 



c 4 y C 

nN , 87rnC s 2 (7 
= 4vr(^ + P r )ovl - + — . (44) 

These two equations represent evolution equations for the quantity — (e + 
47rll) + which represents one of the structure scalars X TF . 



c 4 

The tensor X n R is defined as 



1 



X aP =* R* aips = - Va , ^R* £pps V^V\ (45) 

where R* a ^s = ^VepysRap £p an d Vep~f6 denote the Levi-Civita tensor. Tensor 
X aj3 can be expressed through its trace and trace free part as 

1 1 

X a p = -X T h a p + X TF {x a Xp ~ 2 h <*P)' ( 46 ) 

One can find X TF with the help of the field equations, Eq.(|46]) and ( |2"2"|) -(|211) 
( see for detail [28]) as 

s 2 

X TF = -e - 47rn + — . (47) 
Thus the evolution equations ( H3l and ( l44"l) in terms of Xtf turn out to be 

^TF + 6 c ~ C 5 C A ' 

X TF C , „. A 8nUC s 2 C 
X TF + 3^- = 4vr( M + P r )aA - + — . 
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The corresponding solutions will be 



Xtf = ^j\-^IiC'C 2 + ^ + ^)dr, (48) 

-4tt r f • s 2 
Xtf = -prj [2nC — (fx + P r )ACa]C 2 dt - —. (49) 

The last equation represents that the anisotropy of pressure and electromag- 
netic field affect the initial state of the energy density. 

4.4 Dissipative Geodesic Charged Dust 

Here we consider the case of dissipative geodesic dust to explore the effect 
of dissipation on energy density inhomogeneity. For this purpose, we assume 
that dissipative dust is moving along geodesic, i.e., P r — P±_ — and A—l. 
The corresponding equations take the form 

C ( A 2s 2 \ 12-nqC , , 

(e - 47T/2)' + 3- \e - 4tt/2 + —J - = 0, (50) 

V s 2 C YHqCB ss' , x 

( e _ 47T/1)' = -3-s - 2— ±— + — . (51) 

Equation (J5"T]) leads to a quantity \I/ 

* = £ + / qCBC 2 dr. (52) 



c 3 y 

Equation flBTj) implies that for \1/ = s = 0, we get fx' = similar to the 
uncharged case. For ^ = fx' = 0, we obtain a relation between s and C from 
( 15T|) as discussed in the previous cases. 
Taking fx' = in Eq.( l5T|) . it follows that 

Using Eqs.( )50|) . @ and (!20|) . we obtain the evolution equation of \I/ 

• f2 , „ 47rq' AnqC 2s 2 C 
^-^ = -4^-^ + w + ^, (54) 
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where r 

Q = 12tt [ qCBC 2 dr. (55) 



The general solution of Eq. (I54j) is 



1 , 3 4nq'C 3 , 47rgC 2 C , 2s 2 C , 



^ = c^J Q {n ~ AnafL ~ + g + ~c^ )dt - (56) 

This shows that charged dust density inhomogeneity depends on dissipative 
terms, shear and the electromagnetic field, for the initially homogeneous 
configuration. To investigate the effect of these factors on evolution of we 
consider shearfree case for which Eq.fllj) implies that C = Br. Consequently, 
Eq.flSSD yields 

Q = f ^-dr. (57) 
Jo r 

The corresponding evolution equation is 

* = _L f (Q- Anq'C 2 r + AnqCrC + ^-)dt. (58) 
^ Jo C 

Finally, we discuss the relaxation effects in evolution of The diffusion 
approximation (e = 0) implies q — q, — /x. Thus from Eq.( l2TT) . it follows 
that 

2ss' ,„ n . 

' = -T e + «jci- < 59 » 

Combining the above equation with transport equation (|32|) . we have 



kT 2ss't 



B(l - |0r) AnBC^l - §0r) 
For the shearfree case, this reduces to 

nrT' 2ss'rr 



(60) 



C(l - %Qt) 4ttC 5 (1 - |6r) 



(61) 



When we insert this value of q in Eq.( )58|) . we obtain \l/ in terms of relaxation 
time r which help to analyze relaxation effects in terms of electromagnetic 
field. 
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5 Conclusion 



In this paper, we have studied how charge affects the energy density inhomo- 
geneity and stability of the conformal flatness. For this purpose, the evolution 
equations of the Weyl tensor are formulated. We have investigated different 
aspects of fluid distribution responsible for energy density inhomogeneity. 
The main results are summarized as follows: 

• In the case of non-dissipative charged dust and isotropic fluid, we have 
found that when we take electromagnetic field contribution, the energy 
density inhomogeneity is not controlled by the Weyl tensor alone but it 
also depends on charge. It is worth mentioning here that charge affects 
the conformal flat condition. Here the vanishing of the Weyl tensor 
implies that shear depends on charge. For shearfree fluid, the Weyl 
tensor disappears but this is not true for isotropic fluid. In the absence 
of charge, our result reduce to [24J. 

• For anisotropic fluid, it is found that the shear inhomogeneity also 
depends on charge in addition to anisotropy of pressure. A specific 
quantity associated with energy density inhomogeneity is identified as 
one of the structure scalar. 

• Finally, we have investigated the effect of charge on dissipative geodesic 
dust. The effect of charge and different factors on the evolution of \I/ (a 
quantity that determines the existence of energy density inhomogene- 
ity) are explored. The relaxation effects in the evolution of ^ with the 
inclusion of charge are also indicated in this case. 
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